Population and mass imbalance in atomic Fermi gases. 
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We develop an accurate theory of resonantly interacting Fermi mixtures with both spin and mass 
imbalance. We consider Fermi mixtures with arbitrary mass imbalances, but focus in particular 
on the experimentally available ^Li-^^K mixture. We determine the phase diagram of the mixture 
for different interactions strengths that lie on the BCS side of the Feshbach resonance. We also 
determine the universal phase diagram at unitarity. We find for the mixtures with a sufficiently 
large mass imbalance, that includes the ''Li-*''K mixture, a Lifshitz point in the universal phase 
diagram that signals an instability towards a supersolid phase. 
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I. INTRODUCTION 

Ultracold quantum gases of fermionic atoms arc at 
the center of attention of both experimental and theo- 
retical physicists. Because of the amazing experimen- 
tal control in these gases they offer the possibility to 
experimentally explore various pairing phenomena. As 
a result, many fundamental discoveries have already 
been made, for example the realization of the crossover 
from a Bardeen-Cooper-Schrieffer (BCS) superfluid of 
loosely bound Cooper pairs to a Bose-Einstein conden- 
sate (BEC) of tightly bound molecules, also called the 
BEC-BCS crossover 

For high- Tc superconductors, the size of a Cooper pair 
is comparable to the average distance between the elec- 
trons. This is analogous to the intermediate regime in 
the BEC-BCS crossover. In the BCS regime, the dis- 
tance between two particles making up a Cooper pair is 
much larger than the average distance between particles, 
whereas in the BEC regime the distance between parti- 
cles within a pair is much smaller than the average inter- 
particle distance. If the size of a pair can be manipulated, 
it is possible to go smoothly from one regime to the other. 
This manipulation can be achieved in an atomic Fermi 
gas, where the interaction strength between the atoms in 
the two different spin states can be controlled with a Fes- 
hbach resonance. The smooth BEC-BCS crossover was 
eventuallyrealized in a trapped gas of '^"K [3l| and ^Li 
atoms [Jtl. 

For a Fermi mixture pairing is always possible for an 
equal amount of particles in each spin state. However, 
pairing is absent for the noninteracting system with all 
particles in one spin state. Therefore, as a function of 
population imbalance there must exist a phase transition 
at low temperatures. Experimentally, this transition was 
studied for the strongly interacting and mass-balanced 
case , and the phase diagram was found to be gov- 

erned by a tricritical point that resulted in the observa- 
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tion of phase separation [T^, [HI . Sarma superfluidity is 
also likely to be present in this system Il3i|. but has 
not been unambiguously identified yet. 

Besides an imbalance in the particle densities, there 
can also be an imbalance between the masses of the par- 
ticles. Consequences of these imbalances are interesting 
to study, also because imbalanced Fermi mixtures exist in 
certain condensed-matter systems in a magnetic field, in 
nuclear matter, and even in the quark-gluon plasma that 
is supposed to be present in the core of heavy neutron 
stars 

While studying mass-imbalanced Fermi mixtures, we 
first focused on the mixture with a mass ratio of 6.7, 
corresponding to a ^Li-'*'^K mixture. This resulted in a 
Letter flGl] where we presented the phase diagram of this 
mixture as a function of polarization and temperature in 
the strongly interacting limit, or the so-called unitarity 
limit, where the scattering length of the interparticle in- 
teraction diverges. In the unitarity limit, the size of the 
Cooper pairs is comparable to the average interparticle 
distance and the pairing is a many-body effect. As a 
result, the mass imbalance has profound effects on the 
pairing, because the imbalance strongly affects the two 
Fermi spheres that are present in the system. We showed 
that the phase diagram of the ^Li-'"'K mixture contains 
a Lifshitz point for a majority of heavy fermions 16:]. 
Typically, Lifshitz points are found at weak interactions 
where the critical temperatures are very low. However, 
we found that in the strongly interacting limit the phase 
diagram of the ^Li-'*'^K mixture already contains a Lif- 
shitz point at accessible temperatures, see Fig. [T7] be- 
low. At a Lifshitz point the phase transition undergoes 
a dramatic change of character. Rather than preferring 
a homogeneous order parameter, the system now forms 
an inhomogeneous superfluid. 

The possibility of an inhomogeneous superfluid was 
early investigated by Larkin and Ovchinnikov (LO), who 
considered a superfluid with a single standing-wave or- 
der parameter [l7| , which is energetically more favorable 
than the plane-wave case studied by Fulde and Ferrell 
(FF)[l^. Since the LO phase results in periodic mod- 
ulations of the particle densities, it is a supersolid [19|. 
The FF and LO phases have intrigued the condensed- 
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matter community for many decades, but only very re- 
cently strong evidence for the FFLO phase has been ob- 
tained in a one-dimensional imbalanced Fermi mixture 
of two spin states ^2d\ . Theoretically, it is challenging to 
describe the phase diagram below a Lifshitz point. The 
most stable states are likely to be complicated superpo- 
sitions of standing waves, where different ansatzes lead 
to different stability regions [21, 22]. 

Initially, our focus was on the ^Li-^'^K mixture, be- 
cause this is experimentally a very promising mixture, 
where several accessible Feshbach resonances are identi- 
fied [1^ and both species have also been simultaneously 
cooled into the degenerate regime [l^l . In contrast to our 
previous work, we now consider Fermi mixtures with an 
arbitrary mass ratio. Moreover, apart from the unitarity 
limit, we also cover the BCS regime. We do not calculate 
phase diagrams on the BEC side of the Feshbach reso- 
nance, because in that regime we should include thermal 
molecules in our calculations, which requires a different 
theory. Due to the population and the mass imbalance, 
the two Fermi spheres in the system are typically mis- 
matched, which can induce phase separation with the lo- 
cation of the tricritical point depending on the mass ratio 
and the interaction strength [25[. In this paper, we also 
consider Lifshitz instabilities and also find a multicritical 
point for the mass-balanced case at a very weak interac- 
tion, as shown in Fig. [T31 At this interaction strength 
the Lifshitz point and the tricritical point occur for the 
same polarization and temperature. Furthermore, in the 
imbalanced Fermi mixture the quasiparticle dispersions 
in the superfluid phase also give rise to gapless Sarma 
superfluidity. For the mass-balanced Fermi mixture and 
for the ^Li-'*"K mixture we therefore also calculated the 
regions in the phase diagrams where the superfluid is 
gapless. 

In our Letter [l^ , we showed that although mean-field 
theory vastly overestimates critical temperatures in the 
unitarity limit, it is very useful for a qualitative descrip- 
tion of the physics. From the mass-balanced case, we 
know that the critical temperatures found using mean- 
field techniques are lowered mainly by two effects, namely 
the fermionic selfenergies and the screening of the inter- 
action due to particle- hole fluctuations [2^ [l^]. After 
discussing mean-field theory for the imbalanced Fermi 
gases, we take both these effects into account in the uni- 
tarity limit for the mass-balanced case and for the exper- 
imentally interesting case of the ^Li-''"K mixture. This 
leads to results that compare well with Monte Carlo cal- 
culations [1^ and for equal masses also with experiment 
pj| . Our procedure gives a reduction of the mean-field 
critical temperatures by a factor of 3. This makes it ex- 
perimentally more difficult, but not impossible, to reach 
also for the mass-imbalanced case the superfluid regime. 
Very importantly, the Lifshitz point remains present in 
the phase diagram after taking fluctuations into account. 
This hopefully brings the observation of inhomogeneous 
superfluidity within experimental reach. 

This paper is organized as follows. We start with dis- 



cussing the interactions in Fermi mixtures in Sec. II. In 
Sec. Ill, we give a brief discussion of the Landau theory 
that we use to describe phase transitions. In particu- 
lar, we introduce the Landau thermodynamic potential 
and the order parameter. Next, in Sec. IV, we discuss 
the mean-fleld theory that we use to calculate phase di- 
agrams, which are presented in Sec. IVC. We discuss the 
phase diagrams for three different mass ratios at different 
interaction strengths, to explore the various topologies of 
the phase diagrams that can arise. After that, we discuss 
in more detail the effects of the mass imbalance. We then 
also explore the effect of the interaction strength on the 
position of the tricritical points and Lifshitz points in the 
phase diagram. Subsequently, we discuss the presence of 
the superfluid Sarma phase for the mass-balanced case 
and the ^Li-'^^K mixture, both at unitarity. All these cal- 
culations use in first instance mean-field theory. Then, 
in Sec. V, we include fiuctuation effects to obtain more 
quantitative results for the mass-balanced case and the 
^Li-^°K mixture. We focus here on the unitarity limit, 
although fluctuation effects could also be easily incor- 
porated in the BCS limit of the Feshbach resonance. Fi- 
nally, different appendices are added where more calcula- 
tions can be found. In Appendix A it is explained how the 
thermodynamic potential for the mass-imbalanced Fermi 
gas is obtained, while Appendix B contains the calcula- 
tions for the amplitudes of relevant Feynman diagrams. 

II. INTERACTIONS AND FESHBACH 
RESONANCES 

In this paper we study phase transitions in an imbal- 
anced Fermi gas at different interaction strengths. This is 
in particular relevant if the interaction strength is experi- 
mentally imder control. In atomic Fermi mixtures the in- 
terspecies interaction can be controlled using a Feshbach 
resonance [2^. The effect of the microscopic interaction 
potential can be studied via the two-body transition op- 
erator T^^. The matrix elements of this transition oper- 
ator are directly related to the scattering amplitudes. It 
is defined by 

V\^i+^)^f^^\k), (1) 

where I'tp]^'') the scattering states and V is the mi- 
croscopic interaction potential. To find an expression for 
T^^ we start with the Lippmann-Schwinger equation 

f^^ = V + V K^f^^, (2) 

z- Ho 

where z = E + iO and the notation iO implies the limit ie 
with e J, 0. In this paper we study an imbalanced Fermi 
gas with a point interaction 

l/(x~x') ^FofJlx-x'), (3) 

where Vq is negative, since we are interested in an at- 
tractive interaction. If we now consider the Lippmann- 
Schwinger equation at zero energy ^ = 0, multiply both 
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FIG. 1: The s-wave scattering length a as a function of the 
microscopic interaction strength Vo- Here Vo is negative, since 
we are considering an attractive interaction. 



sides with (k'| from the left and with |k") from the right, 
and we insert a completeness relation in the second term 
on the right-hand side, we obtain 



1 



r2B(o) 



1 



dk 



(27r)3 2e(k)' 



(4) 



where half the reduced kinetic energy, £(k) — ?i^k^/2m, 
is the kinetic energy associated with a mass m that is 
equal to twice the reduced mass, namely 



1 



1 



(5) 



Here, m+ and m_ are the masses of a light and a heavy 
particle, respectively. 

The two-body transition matrix is related to the s-wave 
scattering length a by 



1 



(6) 



Since we are interested in the behavior of the Fermi gas at 
ultralow temperatures, we can use a cut-off momentum 
hK to evaluate the integral in Eq. (|4]) and we obtain, 
using Eq. ©, a relation between the scattering length a 
and the microscopic interaction potential Vb, namely 



rmrVo 



2mAVa + An^h^ ' 



(7) 



This relation is shown in Fig. [TJ It can be seen that for 
small values of Vq the scattering length is negative, which 
means that the Fermi mixture is in the BCS regime. 
Then, for Vb = —{2'iTKf/mK the scattering length di- 
verges, which is called the unitarity limit. For large val- 
ues of Vb the scattering length is positive and the Fermi 
mixture is in the BEC regime. 

For a > the microscopic interaction has a bound 
state with an energy E = —f^jma}. This is the single- 
channel picture of a Feshbach resonance. In this paper 



we focus on the unitarity limit and on the BCS side of the 
resonance. Thus, we look at the case where a diverges 
and at negative scattering lengths a. For those scattering 
lengths we can use the single-channel picture, as long as 
the Feshbach resonance is sufficiently broad. Namely, in 
the limit of a broad resonance the amplitude to be in the 
bare molecular state of the Feshbach resonance turns out 
to be very small fsl, 'so'] . 



III. LANDAU THEORY OF PHASE 
TRANSITIONS 

In order to study the critical behavior of a system, in 
our case an imbalanced Fermi gas, we consider the Lan- 
dau thermodynamic potential density wl(A(x)), with 
A(x) the superfluid order parameter. Near the phase 
transition, where the BCS order parameter A(x) is small, 
the Landau thermodynamic potential density can be ex- 
panded as [3l|, |3^ 



fiAr+.. 



(8) 



where the dots denote the higher orders in |Ap and in 
gradients |VAp. The Landau coefficients in the ther- 
modynamic potential all depend on the temperature and 
on the chemical potentials of the two fermion species. If 
in the thermodynamic potential all coefficients are pos- 
itive, the minimum of the thermodynamic potential is 
located at (A(x)) equal to zero and the system will be 
in the normal state. Whereas a phase transition to a 
superfluid state has occurred when the position of the 
global minimum is located at a nonzero order parameter 
(A(x)), which describes a condensate of bosonic pairs. 
In the case that 7 is positive, it costs energy to have a 
spatially varying superfluid. It is then energetically fa- 
vorable for the system to be homogeneous and therefore 
we can restrict ourselves to a pairing field A independent 
of position. 

We consider first the case where 7 is positive and the 
system is homogeneous. For high temperatures all coef- 
ficients in the thermodynamic potential will be positive 
and the system will be in the normal state. But for low 
temperatures it can occur that certain coefficients change 
sign. Suppose that in the Landau thermodynamic poten- 
tial a is negative and all other coefficients are positive. 
The minimum of the thermodynamic potential will then 
be attained at some nonzero (A) and the Fermi gas will 
be in the superfluid state. Thus, as a changes sign a 
phase transition takes place. The temperature at which 
the transition from the normal state to the superfluid 
state occurs, for given chemical potentials, can be de- 
termined by equating the quadratic coefficient to zero, 
a(Tc) — 0, where is called the critical temperature. 
The phase transition just described is called a second- 
order phase transition and it is characterized by the fact 
that the minimum of the thermodynamic potential shifts 
away from zero continuously, see Fig. [5^. 
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T < Tc 

FIG. 2: The thermodynamic potential density i:i;l(|A|) as a 
function of the order parameter A. Panel (a) shows the behav- 
ior of ljl(1A|) for different temperatures when a second-order 
phase transition occurs and panel (b) when a first-order phase 
transition occurs. 

It is also possible to have a first-order phase transition. 
To explain a first-order phase transition we consider the 
situation where in the thermodynamic potential density 
all coefficients but the fourth-order coefficient /? are pos- 
itive. The thermodynamic potential will then typically 
have two minima. One of these is located at A equal 
to zero and the other one will be located at a nonzero 
value of the order parameter. For higher temperatures 
the minimum located at zero is a global minimum and 
the Fermi gas is in the normal state. If the temperature 
is lowered there will be a point where the two minima 
are equal and for even lower temperatures the minimum 
located at a nonzero order parameter (A) is the global 
minimum. The system is then in the superfluid state. 
The first-order phase transition takes place when these 
two minima are equal, i.e., when ijJl(O) = uj-l{{A)). In 
contrast to a second-order phase transition, the location 
of the global minimum of the thermodynamic potential 
now changes discontinuously from being zero to a nonzero 
(A), see Fig.[2}D. 

Next, we consider the case where 7 is negative. The 
system can then gain energy when the order parameter 
varies in space. Thus, instead of being constant, the order 
parameter will now depend on position. Fulde and Ferrell 
studied the plane-wave solution [l^ 

A(x) = Aoe"'■^ (9) 

while Larkin and Ovchinnikov considered a superfiuid 
with the single standing-wave order parameter [17| 

A(x) = Aocos(k-x), (10) 

which turns out to be energetically more favorable than 
the plane-wave case. Superpositions of more than two 
plane waves are also possible. 

In the LO phase the wavefunction of the bosonic pairs 
is periodic and therefore there will also exist a periodic- 
ity in the atomic density. This periodic structure shows 
itself in the diagonal elements of the one-particle den- 
sity matrix n(x, x') and is therefore called diagonal long- 
range order. This diagonal long-range order is what char- 
acterizes a solid. In a fermionic superfluid a fraction 



of the Cooper pairs is in the lowest energy eigenstate. 
There thus exists a long-range order between the posi- 
tions of the pairs. It is also said that the two-particle 
density matrix g{x, x') has off-diagonal long-range order, 
which implies that ^(x, x') does not vanish in the limit 
|x — x'l — > 00 and characterizes a superfluid. In the case 
where 7 is negative a transition occurs from the normal 
state to a superfluid where the atomic density in the su- 
perfluid has a periodic structure. Then there exists both 
diagonal long-range order as well as off-diagonal long- 
range order and the state is both solid and superfluid. 
This is called a supersolid phase [H, HI] . 

Further on, we show that in a mass-imbalanced Fermi 
gas not only a second-order phase transition can occur, 
but also a first-order transition and even a transition to 
an inhomogeneous superfluid, depending on the values of 
the two chemical potentials. Therefore, there must be 
points in the phase diagram where the character of the 
phase transition changes. If the phase transition changes 
from being second order to first order there will be a 
tricritical point [ssj. This point can thus be found by 
setting 

a(Te3) = /3(Te3) = 0, (11) 

where Tcs is the tricritical temperature. For tempera- 
tures higher than the tricritical temperature the phase 
transition from the normal to the superfluid state is of 
second order, while for T < T^s there is a flrst-order 
phase transition and phase separation occurs. 

The phase transition could also change from being a 
transition from the normal state to a homogeneous su- 
perfluid to a transition from the normal state to a super- 
solid. The point in the pha se diagram where this occurs 
is called a Lifshitz point 33"] . This point can be computed 
by demanding 

a{n) = Jin) = 0, (12) 

where Tl is called the Lifshitz temperature. For temper- 
atures lower than the Lifshitz temperature the transition 
will be from a normal state to a superfluid state where 
the bosonic pairs have a nonzero momentum. Superflu- 
idity at nonzero momentum can be established in many 
ways and due to this variety of possibilities it is difficult 
to predict which kind of superfluidity will be present be- 
low the Lifshitz point. However, they all have to emerge 
from the Lifshitz point and therefore it is important to 
know the position of the Lifshitz point. 

Apart from the above two possibilities, we can think of 
other scenarios for the change in character of the phase 
transitions. But in the phase diagrams we calculated 
for the imbalanced Fermi gas we only found tricritical 
points and Lifshitz points. Therefore, these are the only 
two possibilities we discuss in the following. 



IV. MEAN-FIELD THEORY 

In this section we present the mean-field thermody- 
namic potential, which is an approximation to the exact 
Landau thermodynamic potential for the Fermi gas with 
population and mass imbalance. If we have an expression 
for the Landau thermodynamic potential, we are able to 
determine the phase diagram. 

Although mean-field theory does not contain all in- 
teractions present in a Fermi mixture, it turns out that 
mean-field theory already incorporates all the relevant 
physics determining the topology of the phase diagrams. 
Adding fluctuation effects only c hang es the phase dia- 
grams quantitatively [H, [1^ H^, HJ Because of the 
rather straightforward and transparent calculations in 
mean-field theory, we first discuss mean-field theory in 
some detail. Later on we take fluctuation effects into ac- 
count in order to obtain more quantitative results which 
we can compare with Monte Carlo calculations and for 
the mass-balanced case with experiment. 



A. Thermodynamic Potential 

We consider a two-component Fermi mixture, i.e., a 
mixture containing either a single fermionic species, for 
which two different hyperfine states are present, or con- 
sisting of two different fermionic species with access to 
a single hyperfine state. A balanced Fermi gas consists 
of a single species with an equal population of both spin 
states. In an imbalanced Fermi gas we allow the popula- 
tions and masses to be different. To take into account a 
population imbalance we use different chemical potentials 
for the two (pseudo)spin states, while a mass imbalance 
implies that the particles in the two hyperfine states have 
different masses. The chemical potential and the mass of 
the fermions in state \a) will be denoted by ^„ and mo- 
respectively. From now on heavy particles are always de- 
noted by a minus sign and light particles by a plus sign, 
thus a — ±. 

The mean-field thermodynamic potential for the im- 
balanced Fermi gas is given by 



^l(|A|) 



},(13) 



where the two-body transition matrix r^^(O) is given by 
Eq. ([6]) . This thermodynamic potential is a direct gener- 
alization of the thermodynamic potential for a balanced 
Fermi gas [H, [sl]. In the above expression /3 — l/ksT 
is the inverse thermal energy, /i = -|- /x_)/2 is the 
average chemical potential and half the reduced kinetic 
energy is e(k) = (e+(k) + e_(k))/2 with 
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2ma 




FIG. 3: The single-particle dispersions hwa^y^^ for the bal- 
anced Fermi gas. Panel (a) shows the dispersions for zero A 
where the system is in the normal state. Panel (b) depicts 
the case where the system is in the superfluid state, thus the 
dispersions for nonzero A. In both panels the dashed lines 
are the dispersions of the hole-like excitations, while the full 
lines give the particle-like dispersions. The latter show more 
clearly the opening of a gap 2A at the Fermi level due to the 
formation of a condensate of Cooper pairs. 




k = -Jlm/d/h 



FIG. 4: The dispersions fex)o-(k) for the Fermi gas with both 
population and mass imbalance. Panel (a) shows again the 
dispersions for zero A where the system is in the normal state. 
And panel (b) again depicts the case where the system is in 
the superfluid state, thus the dispersions for nonzero A. The 
meaning of the dashed and full lines are the same as for Fig. 

El 



The first terms in the above thermodynamic potential 
represent the BCS ground state of the mixture, where 
?iw(k) is the average dispersion of the quasiparticles 



(15) 



with |A| the so called BCS gap parameter. The complex 
pairing field A is on average related to the expectation 
value of the pair annihilation operator through 



(A(x)) =Fo(V'-(x)7^+(x) 



(14) 



(16) 



with V'(t(x) the fermionic annihilation operators. 

The second part of the thermodynamic potential, 
namely the part containing the logarithms, corresponds 
to the contribution of an ideal gas of quasiparticles. Here, 
huja{\i) is the dispersion relation of the quasiparticles in 
state |(t), given by 

fiiJo(k) = nu{k) - a[2h - £+(k) + £-(k)]/2, (17) 

with h ~ (fi^ — fJ.-)/2 the difference in chemical poten- 
tials. For the unpolarized Fermi gas with equal masses 
the dispersions fta;o-(k) reduce to the average dispersion 
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(b) 



FIG. 5: (a) The ladder diagram with external momentum q. 
The wiggly lines denote the Cooper pairs, which can break 
up into two fermions of different spin, denoted by the normal 
lines. Here, q is the wavevector of the Cooper pairs, while 
k and q — k are the wavevectors of the fermions. When q 
is equal to zero, the amplitude of this diagram corresponds 
to the quadratic Landau coefficient a in the thermodynamic 
potential, (b) The Feynman diagram corresponding to the 
fourth-order Landau coefficient /?. 



in Eq. This dispersion is plotted in Fig. [3] for 

both the normal state (Fig. [5^) and the superfluid state 
(Fig. [Sjj). In this case the superfluid is gapped and bal- 
anced. For k — |k| < ^2rrni/h the quasiparticle dis- 
persion describes hole-like excitations. If we mirror this 
hole-like part of the quasiparticle dispersion, we obtain 
the negative dispersion of the particle-like excitations. 
For k > ^J2m^/h the quasiparticle dispersion already 
describes the particle-like excitations. 

In the case of a polarized superfluid with equal masses 
the dispersions hujai}^) in Fig. |3] are shifted by the dif- 
ference in chemical potentials 2h. For /i > |A| the dis- 
persion of the majority species becomes negative. When 
this occurs, the occupation of the single-particle states 
associated with the negative part of the quasiparticle 
excitation branch actually lowers the ground-state en- 
ergy. Since this leads to additional majority quasipar- 
ticles and, therefore, additional majority particles and 
minority holes, in the ground state, the ground state be- 
comes a polarized superfluid. The resulting gapless and 
polarized superfluid is called the Sarma phase [sl]. For 
the Fermi gas with both mass and population imbalance 
the dispersions are depicted in Fig. SI where Fig.|3^ shows 
the dispersions for zero A and Fig.|3]D for nonzero A. The 
shape of the dispersions is changed due to the difference 
in mass. It can be seen that in this case also one of the 
dispersions, namely fia;_(k), is negative so that it would 
correspond to a gapless Sarma superfluid. 



B. Landau Coefflcients 

For the imbalanced Fermi mixture, we want to study 
the phase transition from the normal state to the su- 
perfluid state. For that transition we want to obtain a 
phase diagram with the critical temperature as a func- 
tion of the polarization, P = (n+ — nJ)/{n+ + n-) with 
Ucr the density of particles in state \a). 

The particle densities n^. can be determined from the 



thermodynamic potential using [3! 

ac^L(A) 



(18) 



A=(A) 



Moreover, in Sec. lIIII the critical conditions for the transi- 
tion were explained. In mean-field theory, the quadratic 
Landau coefficient a is explicitly given by 



a 



dujL 



1 



9|A|2 
dk 

(2^ V2£(k) 



A|2=0 
1 



T2B(0) 

jV+(k) + jV_(k) - 1 
2(e(k) - 



(19) 



where N^Ck) — l/(exp[/3(ecr(k) — fia)] + 1) are the Fermi 
distribution functions. To determine the temperature of 
the tricritical point, we also need to know the fourth- 
order Landau coefficient /3. It is given by 



13^ 



id\A\^y 



|AP=0 



d^k 



1 



/37V+(k)(7V+(k)-l) 



+ 



1 



e(k) - ^I 



(27r)3 4(£(k) - fiy 
/3iV_(k)(iV_(k)-l) 

(l-iV+(k)-iV_(k)) 



(20) 



Determining 7 from the mean-field thermodynamic po- 
tential is not possible, since we have assumed the bosonic 
pairing field A(x) to be independent of position. Never- 
theless, there is a rather simple way to determine this 
coefficient, using Feynman diagrams 37]. The other co- 
efficients could also have been determined using a dia- 
grammatic language. Namely, the quadratic coefficient a 
corresponds to the so called ladder diagram where the in- 
coming and outgoing bosonic fields A have zero momen- 
tum, see Fig. [5^. Physically, a can be interpreted as be- 
ing proportional to the chemical potential of the Cooper 
pairs. The fourth-order coefficient /3 has a diagrammatic 
representation with four external bosonic fields with zero 
momentum, see Fig. [SJd. 

For the transition to the supersolid phase, we con- 
sider the ladder diagram where the bosonic fields carry 
nonzero momentum q, since the supersolid phase consists 
of bosonic pairs with nonzero momentum. The expres- 
sion for this ladder diagram is 



a(q) 



1 



d^k 



1 



r2B(o) J (27r)3 |2e(k) 

jV+(q-k) + jV_(k)-l | 
e+(q-k)+e_(k)-2^ (■ 



(21) 



The actual shape of a(q) as a function of the exter- 
nal momentum q depends on the values of the differ- 
ent parameters in this expression, such as the chemical 
potentials and the temperature. Depending on those 
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FIG. 6: Depending on the temperature and the chemical potentials, the minimum of a(q) is either at zero or at nonzero 
external momentum hq. In both pictures the mass ratio is r = 6.7. In panel (a), the minimum of Q.{c\) is located at zero 
external momentum. The corresponding polarization is P = —0.61 and the temperature is T /Ty = 0.21 with Tp the Fermi 
temperature. In panel (b), the minimum of Q(q) is attained at a nonzero momentum. Here the polarization is P = —0.7 and 
the temperature is T/T-p = 0.13. 



quantities the minimum of a{c{) is attained either for 
zero or for nonzero external momentum. As explained in 
Sec. lIIIi a second-order phase transition can occur, when 
a quadratic coefficient of the Landau theory changes sign. 
These coefhcients are now given by a{q), where q is the 
wavevector of the bosonic pairs. The sign change occurs 
first for the minimum of a(q), which therefore determines 
whether or not the transition happens at nonzero q. The 
expression for the ladder diagram with nonzero external 
momentum can be expanded in even powers of q 



a(q) = ao + aiq^ + a2q^ 



(22) 



where the dots denote higher order powers in q^. If all 
the coefficients are positive, a(q) has a minimum for 
external momentum zero, so that a transition to the ho- 
mogeneous superfluid phase occurs. But if ai is nega- 
tive, a{q) has a minimum for a nonzero external momen- 
tum and therefore it first becomes zero at some nonzero 
value of q. The minimum of the thermodynamic poten- 
tial will then be located at a nonzero order parameter 
with a nonzero momentum. In other words, it will be 
energetically favorable for the bosonic pairs to have ki- 
netic energy and the phase transition that occurs is a 
transition from the normal state to an inhomogeneous 
superfluid. Comparing this with the Landau theory of 
Sec. Illli we see that in the expansion of a(q), oq can be 
identified with the quadratic coefficient a and ai can be 
identified with 7. Thus, from the ladder diagram with 
external momentum an expression for 7 can be found, 
namely 



9q2 



(23) 



q=0 



Physically, 7 can be interpreted as being proportional to 
the inverse of the effective mass of the Cooper pairs. 



C. Results 

In this section the results using mean-field theory are 
presented. First, we present phase diagrams for three 
Fermi gases with different mass imbalances. Then, we 
study the effect of the mass imbalance on the critical 
temperature for an unpolarized Fermi gas. After this, we 
study the effect of the interaction strength on the tem- 
perature corresponding to a tricritical point or a Lifshitz 
point. Finally, we also consider the superfluid Sarma 
phase. 



1. Phase Diagrams 

With the expressions for the Landau coefRcients, the 
phase diagram can be calculated for a fixed mass ratio 
r and a fixed interaction strength l/kpa. We determine 
the phase diagram as a function of temperature T and 
polarization P = (n+ — n_)/(n4. -|- n_). The mass ratio 
r is given by r = m^/m^. The interaction strength is 
characterized by the s-wave scattering length a and the 
Fermi momentum k-p, which is defined as 



,1/3 



(24) 



where n = ?i+ + n_ is the total particle density. In 
order to obtain a phase diagram independent of the total 
particle density n, we scale the temperature with the 
reduced Fermi temperature 



h^kl 

2m ' 



(25) 



where m is twice the reduced mass, introduced in Eq. ([S]) 
and fee is Boltzmann's constant. 

We present the phase diagram for three different mass 
ratios, namely for r — \, which is the mass-balanced case, 
r = 6.7, which corresponds to a ^Li-^'^K mixtures, and 
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FIG. 7: Phase diagrams of the mass-balanced Fermi gas, r = 1, as a function of temperature T and polarization P at different 
interaction strengths. The temperature is scaled with the Fermi temperature Tp. Panel (a) shows the phase diagram in 
the strongly interacting limit, 1/fcpa = 0. There is a tricritical point (TCP), where the normal state (N), the homogeneous 
superfluid state (SF) and the forbidden region (FR) meet. The shaded area sets the scale for panel (b), where the phase 
diagram for a weaker interaction, 1/fcpa = —1, is shown. Again, there is a tricritical point. There is again a shaded region 
to set the scale for panel (c), but it is too small to see. In panel (c) the phase diagram is shown for a very weak interaction, 
1/fcpa = —3, and the critical temperatures are now extremely low. For this interaction we find a Lifshitz point (LP), below 
which there is an instability towards a supersolid (SS). The size of this supersolid region is not calculated within our theory 
and therefore the dashed lines are only guides to the eye. 




FIG. 8: Phase diagrams of the Fermi gas with mass ratio r = 6.7 for different interaction strengths. Panel (a) shows the phase 
diagram in the strongly interacting limit. The size of the box is the same as for the mass-balanced case such that the effect 
of the mass imbalance can be seen. The shaded region sets the scale for panel (b), where the interaction strength is weaker, 
1/fcpa = —1. And in panel (b) the shaded region in turn sets the scale for panel (c) where the interaction is even weaker, 
1/fcpa = —3/2. The dashed lines are again guides to the eye. 




FIG. 9: Phase diagrams of the Fermi gas with mass ratio r = 10 for different interaction strengths. Again, panel (a) shows the 
phase diagram in the strongly interacting limit and the size of the box is again the same as for the mass-balanced case. The 
shaded region sets the scale for panel (b), where the interaction strength is weaker, 1/fcpa = —1. And in panel (b) the shaded 
region in turn sets the scale for panel (c) where the interaction is even weaker, 1/fcpa = —3/2. Again, the dashed lines are only 
guides to the eye. 
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for r = 10 where an interesting feature in the phase di- 
agram is fomid regarding the tricritical point [2^- For 
these mass ratios we present the phase diagrams at three 
different interaction strengths in order to see what the 
effect of the interaction is on the critical temperature. 
Namely, we present the phase diagrams for a strongly 
interacting Fermi gas with l/fcpa = 0, an intermediate 
interaction strength with 1/k-pa — — 1 and a weakly inter- 
acting Fermi gas. These interaction strengths correspond 
to the BCS side of the Feshbach resonance, see Fig. [T] 

Fig. [Zl shows the phase diagrams of the mass-balanced 
Fermi gas at three different interaction strengths. In 
Fig. [iK the phase diagram for the strongly interacting 
regime where l/fcpa = is depicted. It is symmet- 
ric in the polarizations. The phase transition is of sec- 
ond order for small polarizations. Moreover, we find 
two tricritical points in the phase diagram and below 
the tricritical points the phase transition is of first or- 
der. At a first-order phase transition the order parame- 
ter is discontinuous and therefore also the particle den- 
sities ncr, see Eq. ^S]) . And thus also the polarization 
is discontinuous at a first-order phase transition, which 
gives rise to a forbidden region where phase separation 
occurs. Fig. [7}3 again shows the phase diagram of the 
mass-balanced Fermi gas, but now for a weaker interac- 
tion, namely for l/Zcpa — ~1- Compared to the unitarity 
regime, there are no qualitative changes in the phase dia- 
gram. However, the critical temperatures are lower than 
in the strongly interacting regime. Then, Fig. [Tj: shows 
the phase diagram of the mass-balanced Fermi gas for 
a very weak interaction, namely l/k^a = —3. Apart 
from the fact that the critical temperatures are now ex- 
tremely low, there is a large difference compared to the 
strongly interacting regime, namely we do not find tricrit- 
ical points in the phase diagram but instead we find two 
Lifshitz points. Below the Lifshitz point there is an insta- 
bility towards supersolidity. We assumed that there will 
be a second-order transition from the normal state to the 
supersolid. The critical temperature for this transition is 
found by solving Q;(kLo) = 0; where Iclo is the wavevec- 
tor of the supersolid. The transition from supersolidity to 
the homogeneous superfluid phase is expected to be first 
order [l5l]2]| and thus a forbidden region will be present 
in the phase diagram. This simple scenario is sketched 
in the phase diagram, but it could be that the normal 
to supersolid second-order phase transition is preempted 
by weak first-order transitions from the normal state to 
various more complicated supersolid phases [H, |2l|, [l^j . 
For the calculation of the relevant fourth-order diagrams 
with nonzero external momenta, one needs to make an 
assumption about the crystal structure of the supersolid 
phase, in order to calculate a first-order phase transition 
to the supersolid phase. This calculation for the stabil- 
ity regions of all possible supersolid phases is beyond the 
scope of this paper. To emphasize that we did not in- 
vestigate this region in great detail we there used dashed 
lines in the phase diagram. 

Fig. |5] again shows three phase diagrams at different 



interactions, but now for the mass-imbalanced Fermi gas 
with mass ratio r — 6.7, corresponding to the ^Li-'^^K 
Fermi mixture. In the unitarity regime, see Fig. the 
phase diagram is no longer symmetric in polarizations 
and the temperatures are lower in comparison with the 
mass-balanced case. Furthermore, already in the strongly 
interacting regime we find a Lifshitz point in the phase 
diagram for a majority of heavy atoms. For a majority of 
light particles a tricritical point is found. This is in sharp 
contrast with the mass-balanced case where a Lifshitz 
point is only present for extremely weak interactions. As 
a result, the Lifshitz temperature is then at least about a 
hundred times lower than for the ^Li-'*°K mixture. Below 
the tricritical point there is a forbidden region and below 
the Lifshitz point there is an instability towards a super- 
solid. We sketched the same scenario below the Lifshitz 
point here as for the mass-balanced case. In Fig. [Sja and 
Fig. [S}: the phase diagrams of the Fermi gas are depicted 
at interaction strengths I/Zcfi = —1 and l/fcpa = —3/2 
respectively. The interaction strength does not affect the 
topology of the phase diagram, but it does affect the crit- 
ical temperatures. Just as in the mass-balanced case, the 
critical temperatures are lower for weaker interactions, as 
expected. 

In Fig.|n]we show the phase diagrams of an imbalanced 
Fermi gas with an even larger mass imbalance, namely 
with mass ratio r = 10. The diagram has become even 
more asymmetric. In the strongly interacting regime, we 
also find a Lifshitz point for a majority of heavy atoms. 
But the tricritical point has disappeared, which means 
that the phase transition remains of second order and no 
phase separation occurs for a majority of light particles. 
Below the Lifshitz point the same scenario is sketched 
as before. Fig. shows the phase diagram at a weaker 
interaction strength, namely at l/fcpa = —1. With re- 
spect to the unitarity limit there is a real change in this 
diagram. Namely, the tricritical point reappears in the 
phase diagram. Then, in Fig. I^t the interaction strength 
is even weaker, l/Zcpa = —3/2, and the critical temper- 
atures are lower. But there are no qualitative changes 
with respect to the phase diagram in Fig. ^\p. 



2. Effect of a Mass Imbalance 

By comparing the phase diagrams of the mass- 
balanced Fermi gas and of the Fermi gases with a mass 
imbalance, it can be seen that a mass imbalance causes 
important changes. First, the critical temperatures for 
the mass-imbalanced Fermi gases are lower than for 
the mass-balanced case. Second, for a mass-imbalanced 
Fermi gas the phase diagram is no longer symmetric in 
polarizations. In other words, the maximum critical tem- 
perature is no longer located at zero polarization. Third, 
for Fermi gases with a sufficiently large mass ratio, like 
the two mass-imbalanced Fermi gases we presented here, 
there is a Lifshitz point in the phase diagram for a major- 
ity of heavy atoms. These three changes we now discuss 
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FIG. 10: The critical temperature of an unpolarized Fermi gas 
as a function of the mass ratio r = m-/m+. The circles are 
the numerical results of a = and the full line is the analytic 
result for the weakly interacting limit. Here l/kpa = —3/2. 



FIG. 11: Derivative of the critical temperature with respect 
to the difference in chemical potentials h. The circles are the 
numerical results and the full line is the analytic result. The 
interaction strength is l/fcpa — —3/2. 



in some more detail. 

In the weakly interacting limit, the critical tempera- 
tures are very low. For the unpolarized Fermi gas the 
integral in Eq. can then be evaluated exactly. Then, 
by equating a to zero an analytic result for the critical 
temperature can be obtained |38| 



clp=0 



2v^ Sep 
1 + r Trfce 



Trfcp 



.g7-2g~7r/2feF|a| 



g7-2g-7r/2fcF|a| 



(26) 



where sf is the Fermi energy corresponding to twice 
the reduced mass, fee is Boltzmann's constant and 7 = 
0.5772 is Euler's constant. When the atoms making up 
the Fermi mixture have different masses, r is not equal 
to one. In that case the term (1 + r) is smaller than 
one. And thus, by adding a mass imbalance the criti- 
cal temperature for the unpolarized Fermi gas becomes 
lower. In Fig. [TUl it can be seen that the dependence of 
the critical temperature on the mass ratio is indeed as 
given by Eq. (pS)) and that the agreement with the nu- 
merical results is very good. For the polarized Fermi gas, 
we expect a similar dependence of the critical tempera- 
ture on the mass ratio. 

For the mass-balanced Fermi gas, the highest critical 
temperature is located at zero polarization, see Fig. [71 
whereas for the mass-imbalanced Fermi gas it is located 
at nonzero polarization, see Fig. [5] and Fig. O As a con- 
sequence, the derivative of the critical temperature with 
respect to /i, the difference in chemical potentials, will no 
longer be zero at P = 0. In the weakly interacting limit, 
the coefficient a in Eq. can be expanded around the 
critical temperature Tc and the critical 'Zeeman' field ho 
corresponding to the unpolarized Fermi gas 



aiT, + STc,ho + Sh) ~ 

a{Tc,ho) + -^{Tc,ho)6Tc + ^{Tc,ho)Sh. 



(27) 



For a continuous phase transition, a has to go to zero. 
The term on the left-hand side and the first term on the 
right-hand side are therefore zero. Then we find from the 
other two terms 



STc 
6h 



da/dh 
^ da/dT 
2 r - 1 



= n 



1 



P=0 



1 



1 



^7 — 2g — Tr/2fcF|a 



(28) 



This is indeed nonzero if the mass ratio is not equal to 
one. Eq. is positive for mass ratios larger than one. 
This is in agreement with our findings that the phase 
diagram shifts towards positive polarizations for a mass 
ratio larger than one, see Fig. [5] and Fig. [HI For mass ra- 
tios smaller than one, Eq. (1^51) is negative and the phase 
diagram then shifts towards negative polarizations. In 
both cases the phase diagram shifts towards a majority 
of light particles. Eq. ([^5]) is plotted in Fig. [TTl together 
with the numerical results. The agreement is good, es- 
pecially for small mass ratios. 

For a sufficiently large mass ratio we find a Lifshitz 
point in the phase diagram for a majority of heavy par- 
ticles. To explain why this can be expected we consider 
the kinetic energies of the particles in Eq. ([H)) , which are 
plotted in Fig. [T^ In the mass-balanced case a fermion 
of one species with momentum k has the same energy as 
a fermion with momentum — k of the other species. Cou- 
pling these degenerate states by a condensate of Cooper 
pairs is energetically desirable since the shift of the en- 
ergy levels due to the coupling is now largest, as is well- 
known from the physics of avoided crossings. Therefore, 
it is thus most favorable to form a pair of two fermions 
with the same but opposite momentum, i.e., a fermion 
with momentum k forms a pair with a fermion with mo- 
mentum — k, since then the energy gain as a result of 
pairing is maximal. The pair then has no kinetic energy 
and the pairs form a homogeneous superfluid. 

Typically, the formation of pairs mostly occurs at the 
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FIG. 12: The energies erj{k) of the fermions for a Fermi mix- 
ture with both a population and a mass imbalance, where 
m_ > m+ and n_ > n+. Thus, there is a majority of heavy 
particles. At the Fermi levels, the energy difference between 
particles with different momenta (dashed line) is smaller than 
the energy difference between particles with the same mo- 
menta. 



Fermi energy. There a pair consists of one fermion with 
momentum kp and one with momentum — kp. Since the 
Fermi momentum depends on the density of the particles 
Ua-, see Eq. the Fermi momenta will not be the same 
in the case of a population imbalance. This makes pair- 
ing less ideal than in the unpolarized case. The critical 
temperatures are therefore lower for a polarized mixture, 
see Fig. [71 

In the mass-imbalanced case two fermions with equal 
momentum but from different species do not have the 
same kinetic energy. Thus, if now a fermion of one species 
with momentum k forms a pair with a fermion with mo- 
mentum — k of the other species, there is an energy differ- 
ence that reduces the energy gain that can be obtained 
due to pairing. This explains the fact that the critical 
temperatures are lower for Fermi mixtures with a mass 
imbalance. With both a mass imbalance and a popula- 
tion imbalance pairs can still be formed with two fermions 
with equal but opposite momentum, even though there 
is an energy difference between the Fermi levels. How- 
ever, this occurs only if the mass ratio and the popula- 
tion imbalance are not too large. In the extreme case 
of a large mass imbalance and a large population imbal- 
ance a more favorable scenario is possible, see Fig.[T3 In 
this figure the kinetic energies are plotted for the differ- 
ent fermions, and there is a majority of heavy particles, 
thus A;p__ > fcp^+. It can be seen that the energy dif- 
ference between particles with different momenta is now 
smaller than the energy difference between fermions with 
the same momentum, i.e., it is now energetically more 
favorable to form pairs with a nonzero momentum that 
allows for a direct coupling of the Fermi levels. The first 
point where this is the case is the Lifshitz point. There 



the system becomes a supersolid. 



3. MuHicritical Point 

In the phase diagrams, Fig. [71 Fig. [51 and Fig. [HI we 
find, depending on the mass ratio and the interaction 
strength, tricritical points (PcSjTcs), which are solutions 
of Eq. ([TT|) . as well as Lifshitz points (PljTl), solutions 
of Eq. (HH). 

At unitarity, we find for the mass-balanced case at uni- 
tarity a tricritical point, whereas we find a Lifshitz point 
for the ^Li-^^K mixture for a majority of heavy parti- 
cles. This Lifshitz point is a clearly distinct point from a 
tricritical point, just as the tricritical point for the mass- 
balanced case is a clearly distinct point from a Lifshitz 
point. This is in accordance with the findings of Parish et 
al. [HI, who found that for a mass-balanced Fermi mix- 
ture the FFLO line detaches from the tricritical point 
away from the BCS limit. However, there is a mass ratio 
between r — I and r — 6.7 for which the phase diagram 
at unitarity contains a point which is both a tricritical 
point and a Lifshitz point. In other words, for this mass 
ratio Eq. (fTTj) and Eq. (|12p have a solution exactly at the 
same point in the phase diagram. This is a multicritical 
point and it is found to occur for the mass ratio r = 4.22. 

It can also be seen in the different phase diagrams that 
for a given mass ratio the location of tricritical and Lif- 
shitz points changes as the interaction strength changes. 
In Fig. [131 it is shown for the mass-balanced Fermi gas 
and for the Fermi gas with mass ratio r = 6.7 how the 
temperature T of these points changes as a function of 
the interaction strength. Note that the polarization is 
not constant for the lines in Fig. [T^l In both Fig. [T5k 
and Fig. 113b one can recognize an exponential decay as 
a function of the interaction strength, which is exactly 
what one would expect from Eq. ([^5|. 

Fig.[T5b shows the mass-balanced case. In the unitarity 
regime and for small values of l/fcp|a| a tricritical point 
is present in the phase diagram, with temperature Tcs- 
This is the full line in Fig.[T5k. Then for some value of the 
interaction strength Eq. PT|) and Eq. have a solution 
at the same point. Thus, there is a multicritical point and 
the temperature of this point is in Fig. [T^ denoted by Tm- 
For even weaker interactions it then turns out that there 
is an instability towards a supersolid and thus we find 
a Lifshitz point in the phase diagram, with temperature 
Tl. This is the dashed line in Fig. [T3k . 

Although the tricritical and Lifshitz points are very 
close together for very weak interactions, in our ap- 
proach, where we perform the momentum integrals to 
calculate the coefficients 7 in Eq. and (3 in Eq. (|20p , 
they are essentially always two distinct points. Except 
for one value of the interaction strength where we find 
a multicritical point. If one, however, assumes particle- 
hole symmetry in the weakly-interacting limit, Eq. (j23p 
and Eq. (1201) become, up to a constant, the same equa- 
tion and thus one finds that the tricritical point and the 
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FIG. 13: Temperature of the tricritical points and the Lifshitz points as a function of the interaction strength. Panel (a) shows 
the mass-balanced case. In panel (b) we look at a Fermi gas with mass ratio r — 6.7 and a majority of heavy particles. Tcs is 
the temperature of a tricritical point, Tm is the temperature of a multicritical point and Tl of a Lifshitz point. 



Lifshitz point are the same point [l^ . 

Fig. [T?b gives the temperature of the Lifshitz point 
present in the phase diagram of the Fermi gas with mass 
ratio r = 6.7 for a majority of heavy particles. Here, we 
do not find a change in character of the phase transition 
as a function of interaction strength. The Lifshitz point 
remains a Lifshitz point. For a majority of light particles 
of this Fermi mixture, we found a tricritical point in the 
strongly interacting limit. We expect that this point will 
also change to a Lifshitz point for very weak interactions, 
but the extremely low temperatures make it numerically 
difficult to investigate this possibility in detail. 



4- Sarma Phase 

In Sec. IVA, while explaining the dispersions of the 
quasiparticles, we discussed the possibility of a Sarma 
phase, which occurs when one of the dispersions hujaO^) 
becomes negative. At each point in the phase diagram 
the value of the order parameter | A| can be found by min- 
imizing the thermodynamic potential in Eq. (jl3p . With 
this value of the order parameter it can be determined 
whether the dispersions are always positive or become 
negative for some range of momenta. In the superfluid 
phase it turns out that there are different regions. In 
one region the dispersions are always positive and the 
superfluid is gapped, while in the other region one of the 
dispersions becomes negative and we are dealing with 
a gapless superfluid. At nonzero temperatures, there 
is no phase transition between these regions. Rather, 
there is only a smooth crossover. We calculated the 
gapless superfluid region in the unitarity limit for the 
mass-balanced case and for the ^Li-^°K Fermi mixture. 
The results are shown in Fig. [Mj In the mass-balanced 
case, Fig. [Hk . there are two regions in the phase dia- 
gram where we have a gapless superfluid. In Fig. 114b the 
mass-imbalanced case is shown. Here, the Sarma region 



is not symmetric in polarizations, just as the rest of the 
phase diagram. The Sarma region is very large for pos- 
itive polarizations, whereas for negative polarizations it 
becomes very small. 

Fig.lTSjdcpicts the distribution functions iVj'^(k) of the 
heavy and the light particles for some point in the phase 
diagram of FiglTib that lies in the superfluid Sarma re- 
gion. In the superfluid phase the distribution functions 
are modified with respect to the Fermi distribution func- 
tions A'^cr(k) and are given by 

It can be seen in Fig. [TSl that the distribution function of 
the heavy particles is nonmonotic, which is a signature 
of the Sarma phase. 

V. FLUCTUATION EFFECTS 

As argued before, mean-field theory only leads to a 
qualitative description of the phase transitions that oc- 
cur in an imbalanced Fermi gas with unitarity-limited 
interactions. To achieve also a quantitative description, 
we have to take fluctuation effects into account. From 
renormalization-group calculations [13] , we know that es- 
pecially selfenergy effects and the screening of the inter- 
action by particle-hole fluctuations are important in the 
unitarity limit. The corresponding corrections result in 
quantitative agreement with experiments in the mass- 
balanced case [27| , and give rise to more accurate predic- 
tions for upcoming experiments with the very promising 
^U-^°K mixture [l3. 

As explained in Refs. [l3,[l^, a convenient way to take 
selfenergy corrections into account is to introduce renor- 
malized chemical potentials fj,'^ that describe the selfen- 
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FIG. 14: Phase diagrams at unitarity including the region where the superfluid is gapless. Panel (a) shows the mass-balanced 
case. It is the same phase diagram as in Fig. 8a, only now the superfluid region is divided into a gapped superfluid region 
(BCS) and two Sarma regions (S). Panel (b) shows the mass-imbalanced Fermi gas with mass ratio r = 6.7. The small cross in 
panel (b) is the point at which we calculated the distribution functions of the particles shown in Fig. 1151 




FIG. 15: The distribution functions A^f (k) for the heavy 
particles, panel (a), and the light particles, panel (b), in the 
gapless superfluid phase. 



ergy of particles with spin cr in a Fermi sea of particles 
with spin —a. This can be achieved by using 



/4 + M-a ■ 



(30) 



where Cc is a coefficient that can be determined in 
the ladder approximation [ill, but also with the use of 
renormalization-group calculations 27[ or Monte Carlo 
calculations [31. For the mass- balanced case, we use 
c± = 0.6, while for the ^Li-'^'^K mixture we use c+ = 2.3 
and c_ = 0.36 to incorporate the Monte Carlo results. 
The substitution of these renormalized chemical poten- 
tials fi'^l^a-) in the mean- field thermodynamic potential 
uji,{A; ^'^) results in the following equation for the den- 
sities 



dfia 



(31) 



A = (A> 



In Ref. [l6|, an accurate comparison between the use of 
renormalized chemical potentials and Monte Carlo calcu- 



lations was made at zero temperatures, leading to excel- 
lent agreement. Note that in the superfluid state Eq. ([50)) 
overestimates the effects of the interactions on the renor- 
malization of the chemical potentials and another cor- 
rection proportional to has to be subtracted in the 
right-hand side of Eq. §UI 

Fluctuations do not only affect the selfenergies of the 
fcrmions in the normal state. There is another effect of 
particle-hole fluctuations that affects the transition to the 
superfluid state. Namely, there is a change in the coef- 
ficient a due to screening of the interspecies interaction, 
which is also called the Gor'kov correction. We can take 
the screening into account by considering an effective 
two-body interaction that includes the so-called random- 
phase approximation (RPA) bubble sum. This procedure 
is diagrammatically represented in Fig. 1161 The inclusion 
of the infinite geometric series of bubble diagrams leads 
at zero external momentum and frequency to 



1 

vZ 



1 



m(o,o), 



(32) 



where hll{0,0) is the amplitude of the bubble diagram. 
In Appendix B, it is shown that this amplitude is given 
by 



m(o,o) 



(27r)3 2/i'-e+(k') + e-(k') 



(33) 



where we use renormalized chemical potentials to also in- 
clude the fermionic selfenergy effects. Using the screened 
interaction of Eq. (|32l) in Eq. Q , we obtain a screened 
two-body transition matrix T^^{0), which consequently 
enters the expression for the quadratic coefficient in Eq. 
p9t . When the quadratic coefficient including screen- 
ing becomes zero, a second-order transition can occur, so 
that the critical condition now becomes 



(34) 
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FIG. 16: The screened interaction containing an infinite sum of bubble diagrams. Panel (a) shows the perturbative expansion 
of the screened interaction and panel (b) shows the result of the resummation of this expansion. 
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FIG. 17: Phase diagrams of strongly interacting Fermi gases, i.e., l/fcpa = 0, where fluctuation effects have been taken into 
account. Panel (a) shows the mass-balanced case, i.e., r = \. The open squares are data along the phase boundaries from the 
experiment of Shin et al. [ll||. The dashed lines denote the first-order phase transitions, not calculated within our theory and 
therefore only guides to the eye. However, they include the correct critical polarization at zero temperature, denoted by a cross, 
known from Monte Carlo calculations [ij]. To calculate the first-order phase transition lines more accurately with fluctuation 
effects taken into account other methods are needed ^1^. Panel (b) shows the phase diagram for the ''Li-*''K mixture, where 
r = 6.7. For both mixtures the critical temperatures are lowered by taking fiuctuation effects into account, compare Fig. [7^ 
and Fig. [Dj,. 



The new critical condition that inchides both screening 
and fermionic selfenergy effects typicaUy reduces the ob- 
tained critical temperatures with a factor of three. 

If we apply this procedure to determine the hne of 
second-order phase transitions in the mass-balanced case, 
we obtain the result in Fig. [T7k . In this figure, also 
the data along the phase boundaries from the experi- 
ment of Shin et al. [ll| are shown. For the unpolarized 
Fermi gas, we find = O.lSTp and ^(Tc) = O.SlTp 
[l^ . which is to be compared with the Monte Carlo re- 
sults Tc = O.lSTp and /i(Tc) = 0.49rp ;2S]. Moreover, 
for the location of the mass-balanced tricritical point we 
find fceTca = 0.09ep+ and Pc3 = 0.24 which is rather 
close to the experimental data [llj . We included the crit- 
ical polarization at zero temperature known from Monte 
Carlo calculations [i^, in order to be able to sketch the 
first-order line and the forbidden region. Other methods 
are needed to calculate the first-order phase transition 
more accurately with fluctuation effects taken into ac- 
count 13]. In Fig. [T7b . we show the line of second-order 
phase transitions for the ^Li-'*^K mixture, when fiuctu- 
ation effects are included. Compared to the mean-field 
result of Fig. [Hb . the critical temperatures are signif- 
icantly lower. For the positions of the tricritical point 
and the Lifshitz point, we use Eq. (PU)) and Eq. (1^ . 
where we again insert the renormalized chemical poten- 
tials to include the fermionic selfenergy effects. We then 



find /cbTcs = O.OSTp and Pc3 = 0.47, and fceTL = 0.05Tp 
and Pl — —0.18, respectively [l^. It is important to 
note that although the fluctuations have quantitatively 
a very large effect, the topology of the phase diagrams 
remains the same. 



VI. DISCUSSION AND CONCLUSION 

In this paper we considered Fermi mixtures consist- 
ing of two different species of fermions, which can both 
have a population and a mass imbalance. On the mean- 
fleld level we calculated the phase diagrams for those 
Fermi mixtures as a function of temperature and polar- 
ization. We calculated the phase diagrams for different 
mass ratios and for different interaction strengths, where 
we found that a mass imbalance leads to a phase diagram 
that is asymmetrical in the polarization. We also consid- 
ered the possibility of a Lifshitz instability. We found 
such instabilities in the Fermi mixtures with a mass ra- 
tion of r = 6.7 and r = 10 in the unitarity limit. For a 
mass-balanced mixture, Lifshitz points occur only in the 
weakly interacting regime. 

By studying the effects of a mass imbalance in more de- 
tail we found analytic results for the critical temperature 
and the change in critical temperature at zero polariza- 
tion. Furthermore, we investigated what happens to the 
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position of the Lifshitz points and the tricritical points 
when changing the interaction strength. In the mass- 
balanced case, we also found a multicritical point for 
weak interactions. Both for the mass-balanced and the 
^Li-^'^K Fermi mixtures, we calculated the regions where 
the superfiuid phase is gapless in the unitarity limit, i.e., 
where the mixtures are in the Sarma phase. These re- 
gions were present at nonzero temperatures, and turned 
out to be quite large. 

Finally, to obtain more quantitative results we intro- 
duced renormalized chemical potentials that include self- 
energy effects to account for the resonant interactions. 
We also took screening effects on the critical tempera- 
ture into account. In this way, we obtained a phase dia- 
gram for the mass-balanced mixture that agrees well with 
Monte Carlo calculations and with experiment. And we 
hope to have obtained a good quantitative description of 
the phase diagram for the ^Li-^"K mixture, where espe- 
cially the presence of a Lifshitz point is exciting. Below 
the Lifshitz point various supersolid states are compet- 
itive. This leads to open questions for further research 
on the interesting possibility of an inhomogeneous super- 
fluid. 
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Appendix A: Derivation of the Mean-Field 
Thermodynamic Potential 

In this Appendix we explicitly derive the mean-field 
thermodynamic potential in Eq. We start with the 

microscopic action for an interacting Fermi mixture con- 
sisting of fermions present in two hyperfine states. It is 
given by 



h/3 







dx 



+ Fo<^;(x, t)4>*_ (x, t)0_ (x, r)0+(x, r)} , (Al) 

where Vq (5(x — x') is the bare point-like interaction asso- 
ciated with the short-range atomic potentials. The last 
term in this action is a fourth-order term in the fermionic 
fields. Therefore, the integral over the fermionic fields in 
the partition function 



cannot be performed analytically. To deal with 
this fourth-order term we introduce by means of a 
Hubbard-Stratonovich transformation bosonic pairing 
fields A(x, r), which are on average related to the 
fermionic fields 0cr(x, r) as in Eq. p6)) . The Hubbard- 
Stratonovich transformation is performed by inserting 
into the partition function the following identity 



d[A*]d[A]e('^"'^+"^-^«l^o '|A*-Vo010;)/ft^ (A3) 



where the inner product in the exponent is a short hand 
notation for 

At j dx(A*(x,T)-0;(x,T)C(x,T)Fo) 

V'(A(x,t) - Fo0_(x,r)0+(x,r)). (A4) 



Inserting Eq. (jA3| into Eq. (jA2|) leaves us with an ac- 
tion S'[A*, A, 0*, (/)] which depends only quadratically on 
the fermionic fields (/)cr(x, r). It is given by 



5[A*,A,, 



i 1 I'd 



-hV] / dr' / dx'C(x,T)GoTi(x,T;x',r')0.(x',T') 

-I- r+ (x, T)r- (x, r) A(x, r) + A* (x, r)</)_ (x, t)0+ (x, t) | . 

(A5) 

The noninteracting Green's function is given by 
G„^i(x,r;x',r') 
= -|.--^-,.|.(x-x>(.-A (A6) 

Note that by performing the Hubbard-Stratonovich 
transformation we have introduced the order parameter 
A(x, r) in an exact manner into the many-body theory. 
If we now integrate out the fermionic fields we will be 
left with an effective action 5°*^ [A*, A], which is directly 
related to the Landau free energy. If we apply mean-field 
theory we assume that the bosonic fields are position and 
time independent, i.e., A(x, r) = A. In order to integrate 
out the fermionic fields, we need to diagonalize the ac- 
tion. First, Fourier transforming the above action and 
then writing it in a more compact way using matrix mul- 
tiplication yields 



Vo 



hj2 [<^:,+(k),0_„,_(-k)] G^i(k,ic.„ 



Z= / d[r]d[0]exp{-5[0*,<^]/n} 



(A2) 



0„.+ (k) 
'-„,-(-k 

(A7) 
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where V is the volume. In the above expression 

a + e+{k) - fi+ 

A* -{ihijj. 



we find to be 



-ihuj„ 



A 

£_(k)-/i_) 



(A8) 



where w„ are the odd Matsubara frequencies. By writ- 
ing the action in matrix form, we have interchanged the 
fermionic fields 4>- and (^1 and thereby we have picked up 
a constant term, namely the sum I]k(e-(k) — ^_) in the 
first line of Eq. (jA7l) . We can now diagonalize this action 
by means of a Bogoliubov transformation. This trans- 
formation consists of unitarily transforming the atomic 
fields (j)n.,aO^) to the quasiparticle fields ■i/'ri,cr(k) 



Vn,+ (k) 




>-n,-(-k). 





u(k) -i;(k) 
v*{k) u*{k) 



where this transformation is unitary if 

Kk)p-f Kk)p = i. 



-„,-(-k) 



(A9) 



(AlO) 



The coefficients of this transformation can be found by 
demanding the off-diagonal matrix elements to be zero. 
The action in terms of the fields ipn,aO^) then reads 



Vn 



hl3}_^{e{k)-fi-hLo{k)) 

k 

■ ^ [{-ihwn + fu^+{k))tp^ _^_(k}■^|Jn^+{k) 



k.n 



+ (-ifii^„ + ricj_(k))V':^_(k)V'„^_(k)] 



(All) 



where the extra term fta;+(k) inside the first sum 
again comes from interchanging fermionic fields. The 
mean-field partition function reads after performing the 
Hubbard-Stratonovich transformation and the Bogoli- 
ubov transformation 

= exp {PV^-^ + Tr[log(-G-i)]} 



= exp{-S''^[A*,A]/h}. 
where G^^ is 

ihuJn — /lW-|_(k) 



hG- 







ihuJn — hid- (k) 



(A12) 



(A13) 



The interaction Vq in the effective action can be elim- 
inated in favor of the two-body transition matrix using 
Eq. Q . It is straightforward to calculate the mean-field 
thermodynamic potential from the effective action, which 



coL{\A\) = -—\ogZ = 



T2B(0) 

lAP 



ij2^og[l + e~^'^-'(^^)}, (A14) 



2e(k) 

where to obtain this thermodynamic potential the sum 
over the Matsubara frequencies was performed using the 
following identity 

log(/3(-ifii^„ + e)) = log (1 + e-^') , (A15) 

n 

which can be derived using contour integration. 

Appendix B: Feynman Diagrams 

To determine the critical temperature at which a phase 
transition will take place we can, as mentioned in Section 
II, either calculate coefficients from the thermodynamic 
potential by deriving it with respect to |A| or we can 
use Feynman diagrams. We present the latter scheme in 
some detail here. 

From the action in Eq. (|A5|) we can determine the 
propagator for the fermionic fields and the vertices for 
the interactions between the fermions and the bosons. 
Namely the propagator in momentum space is 

Go,CT(k,iw„) = — , / (Bl) 

This propagator is in the Feynman diagrams in Figs. [SJi 
and [5)3 represented by a straight line, where the plus 
(minus) indicates a light (heavy) particle. The inter- 
action vertex is proportional to 5(k± + k2 + k3)Sn,n' ,n" , 
where k^ are the momenta and n determine the frequen- 
cies of the incoming and outgoing particles. This repre- 
sents nothing but conservation of momentum and energy. 
In the diagrams in Figs. [S^ and [5}d the vertices are the 
points where three propagators meet. 

In order to determine the Lifshitz point we need an 
expression for the ladder diagram with nonzero external 
momenta. 

1 f dk 

_ 1 ^ /■ dk -h 

'm^J (27r)3iw„ft-H(£+(q-k)-/i+) 

-h 



-iujnh + (e-(k) - ^_) 



(B2) 



We can split the fractions and then perform the summa- 
tion over the Matsubara frequencies. This results in 

dk N+{ci-k) + N^{k)^l 



(27r)3 £+(q-k) + e_(k)-2/i' 
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If we want to obtain the full expression for the quadratic 
part, we have to add the terms proportional to |Ap 
from the effective action. By doing so we obtain a{q) 
in Eq. (I2ip . In the same way the diagrams can be calcu- 
lated needed for a in Eq. ([T9)) and for (3 in Eq. ([20)) . 

When we include fluctuation effects, we also take the 
screening of the interaction by the particle-hole fluctua- 
tions into account. We do so by replacing the bare in- 
teraction potential Vq by a screened interaction potential 
Vsc containing the infinite sum of so-called RPA bubble 
diagrams, see Fig. [Tni The sum over all the bubble dia- 
grams in Fig. ll6l is a geometric series. Using this, we find 
the result for the screened interaction Eq. (|32|). Thus, 
to find the screened interaction, we have to calculate the 
amplitude of the bubble diagram. The bubble diagram 
consists of two fermionic propagators with momentum 



and energy going in opposite directions through the dia- 
gram. The amplitude is given by 



/■_dk 



— ifujJn + £-(k) — 



(B3) 



Writing the above expression as the sum of two fractions 
and then performing the sum over the Matsubara fre- 
quencies results in 

^ , , /■ dk iV+(k) - A^_(k) , , 
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